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^y ics of semirings to numerical analysis and computing. Concepts 

of universal algorithm and generic program are discussed. Rela- 



tions between these concepts and mathematics of semirings are 
examined. A very brief introduction to mathematics of semir- 
ings (including idempotent and tropical mathematics) is presented. 
Concrete applications to optimization problems, idempotent linear 
algebra and interval analysis are indicated. It is known that some 
nonlinear problems (and especially optimization problems) become 
linear over appropriate semirings with idempotent addition (the 
so-called idempotent superposition principle). This linearity over 
semirings is convenient for parallel computations. 
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1. Introduction 

It is well known that programmers are the laziest persons in the 
world. They constantly try to write less and get more. The so-called art 
of programming is just the way to achieve this goal. They started from 
programming in computer codes, found this practice boring and in- 
vented assembler, then macro assembler, then programming languages 
to reduce the amount of work (but not the salary) . The next step was to 
separate algorithms and data. The new principle "Algorithm + data 
structure = program" was a great step in the same direction. Now 
some people could work on algorithms and express those in a more 
or less data structure independent manner; others programmers im- 
plement algorithms for different target data structures using different 
languages. 

This scheme worked, and worked great, but had an important draw- 
back: for the same algorithm one has to write a new program every 
time a new data type is required. Not only is this a boring process, 
it's also consuming time and money and worse, it is error-prone. So 
the next step was to find how to write programs in a way indepen- 
dent of data structures. Object oriented languages, such as C++, Java 
and many others (see, e.g., [47, 62]) opened precisely such a way. For 



MATHEMATICS OF SEMIRINGS 3 

C++, templates and STL give even more opportunities. This means 
that for a given algorithm one can vary a data structure while keeping 
the program unchanged. The principle was changed to "Algorithm + 
data structure family = program." 

But how does this approach work? What constitutes a "data struc- 
ture family"? The answer is this: operations. Any algorithm manip- 
ulates data by means of a set of basic operations. For example for 
sorting the comparison operation is required; for scalar products and 
matrix multiplications, the operations of addition and multiplication 
are required. So to use a new data structure with the same "generic" 
program one has to implement the required operations for this data. 

The hope was that some standard libraries of generic programs would 
cover almost all programmers' needs and programming would be very 
much like Lego constructing. Well, to some extent it worked, but not 
as much as it was expected. Why? For different reasons. Some of them 
are economical. Why invest time and money in solving generic needs 
when a fast patch exists? Indeed, who would buy the next release 
if the current one is perfect? Another reason: there is not a great 
variety of possible data structures for most popular algorithms. In 
linear algebra one can use floating numbers, double precision numbers, 
infinite precision floating numbers, rational numbers, rational numbers 
with fix precision, Hensel codes, complex numbers, integers. Not much. 

Mathematics of semirings gives a new approach to the generic pro- 
gramming. It parameterized algorithms. The new principle is "Pa- 
rameterized algorithm + data structure family = program." What are 
these parameters? They are operations (e.g., addition and multiplica- 
tion). Sounds great, but does it really work? Yes. And we will show 
how and why. For example, we will show how the same old algorithm 
of linear algebra transformed by a simple redefinition of operations, 
can be applied to different problems in different areas. For example, 
the same program for solving systems of linear equations can be ap- 
plied to the shortest path problem and other optimization problems 
(and interval versions of the problems) by means of simple redefini- 
tions of the addition and multiplication operations. The algorithm was 
changed (by changing parameters, i.e. operations), the data structure 
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was changed, the nature of the problem was changed, but the program 
was not changed! 

There are deep mathematical reasons (related to mathematics of 
semirings, idempotent and tropical mathematics) why this approach 
works. We will briefly discuss them in this paper. 

The concept of a generic program was introduced by many authors; 
for example, in [36] such programs were called 'program schemes.' In 
this paper, we discuss universal algorithms implemented in the form 
of generic programs and their specific features. This paper is closely 
related to papers [38-40, 42, 43], in which the concept of a universal al- 
gorithm was defined and software and hardware implementation of such 
algorithms was discussed in connection with problems of idempotent 
mathematics [31, 32, 37-46, 48-53, 55, 71, 72]. In the present paper 
the emphasis is placed on software and hardware implementations of 
universal algorithms, computation with arbitrary accuracy, universal 
algorithms of linear algebra over semirings, and their implementations. 

We also present a very brief introduction to mathematics of semir- 
ings and especially to the mathematics of idempotent semirings (i.e. 
semirings with idempotent addition). Mathematics over idempotent 
semirings is called idempotent mathematics. The so-called idempotent 
correspondence principle and idempotent superposition principle (see 
[38-43, 49-53] are discussed. 

There exists a correspondence between interesting, useful, and im- 
portant constructions and results concerning the field of real (or com- 
plex) numbers and similar constructions dealing with various idempo- 
tent semirings. This correspondence can be formulated in the spirit 
of the well-known N. Bohr's correspondence principle in quantum me- 
chanics; in fact, the two principles are intimately connected (see [38- 
40]). In a sense, the traditional mathematics over numerical fields can 
be treated as a 'quantum' theory, whereas the idempotent mathematics 
can be treated as a 'classical' shadow (or counterpart) of the traditional 
one. It is important that the idempotent correspondence principle is 
valid for algorithms, computer programs and hardware units. 
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In quantum mechanics the superposition principle means that the 
Schrodinger equation (which is basic for the theory) is linear. Simi- 
larly in idempotent mathematics the (idempotent) superposition prin- 
ciple (formulated by V. P. Maslov) means that some important and 
basic problems and equations that are nonlinear in the usual sense 
(e.g., the Hamilton- Jacobi equation, which is basic for classical me- 
chanics and appears in many optimization problems, or the Bellman 
equation and its versions and generalizations) can be treated as linear 
over appropriate idempotent semirings, see [38-41, 48-51]. 

Note that numerical algorithms for infinite-dimensional linear prob- 
lems over idempotent semirings (e.g., idempotent integration, integral 
operators and transformations, the Hamilton-Jacobi and generalized 
Bellman equations) deal with the corresponding finite-dimensional ap- 
proximations. Thus idempotent linear algebra is the basis of the idem- 
potent numerical analysis and, in particular, the discrete optimization 
theory. 

B. A. Carre [10, 11] (see also [7, 23-25]) used the idempotent linear 
algebra to show that different optimization problems for finite graphs 
can be formulated in a unified manner and reduced to solving Bellman 
equations, i.e., systems of linear algebraic equations over idempotent 
semirings. He also generalized principal algorithms of computational 
linear algebra to the idempotent case and showed that some of these 
coincide with algorithms independently developed for solution of op- 
timization problems; for example, Bellman's method of solving the 
shortest path problem corresponds to a version of Jacobi's method for 
solving a system of linear equations, whereas Ford's algorithm cor- 
responds to a version of Gauss-Seidel's method. We briefly discuss 
Bellman equations and the corresponding optimization problems on 
graphs. 

We stress that these well-known results can be interpreted as a man- 
ifestation of the idempotent superposition principle. 

We also briefly discuss interval analysis over idempotent and positive 
semirings. Idempotent internal analysis appears in [45, 46, 69]. Later 
many authors dealt with this subject, see, e.g. [12, 20, 27, 57, 58, 77]. 
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It is important to observe that intervals over an idempotent semiring 
form a new idempotent semiring. Hence universal algorithms can be 
applied to elements of this new semiring and generate interval versions 
of the initial algorithms. 

Note finally that idempotent mathematics is remarkably simpler 
than its traditional analog. 

2. Universal algorithms 

Computational algorithms are constructed on the basis of certain 
primitive operations. These operations manipulate data that describe 
"numbers." These "numbers" are elements of a "numerical domain," 
i.e., a mathematical object such as the field of real numbers, the ring of 
integers, or an idempotent semiring of numbers (idempotent semirings 
and their role in idempotent mathematics are discussed in [10, 11, 22- 
26, 31, 32, 37-43] and below in this paper). 

In practice elements of the numerical domains are replaced by their 
computer representations, i.e., by elements of certain finite models of 
these domains. Examples of models that can be conveniently used 
for computer representation of real numbers are provided by various 
modifications of floating point arithmetics, approximate arithmetics of 
rational numbers [44], and interval arithmetics. The difference between 
mathematical objects ("ideal" numbers) and their finite models (com- 
puter representations) results in computational (e.g., rounding) errors. 

An algorithm is called universal if it is independent of a particular 
numerical domain and/or its computer representation. A typical ex- 
ample of a universal algorithm is the computation of the scalar prod- 
uct (x, y) of two vectors x = (xi, . . . , x n ) and y = (yi, . . . , y n ) by the 
formula (x,y) = X\y\ + • • • + x n y n . This algorithm (formula) is in- 
dependent of a particular domain and its computer implementation, 
since the formula is well-defined for any semiring. It is clear that one 
algorithm can be more universal than another. For example, the sim- 
plest Newton-Cotes formula, the rectangular rule, provides the most 
universal algorithm for numerical integration; indeed, this formula is 
valid even for idempotent integration (over any idempotent semiring, 
see below and [4, 32, 38-43, 48-51]. Other quadrature formulas (e.g., 
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combined trapezoid rule or the Simpson formula) are independent of 
computer arithmetics and can be used (e.g., in the iterative form) for 
computations with arbitrary accuracy. In contrast, algorithms based 
on Gauss-Jacobi formulas are designed for fixed accuracy computa- 
tions: they include constants (coefficients and nodes of these formulas) 
defined with fixed accuracy. Certainly, algorithms of this type can 
be made more universal by including procedures for computing the 
constants; however, this results in an unjustified complication of the 
algorithms. 

Computer algebra algorithms used in such systems as Mathematica, 
Maple, REDUCE, and others are highly universal. Most of the stan- 
dard algorithms used in linear algebra can be rewritten in such a way 
that they will be valid over any field and complete idempotent semiring 
(including semirings of intervals; see below and [45, 46, 69], where an 
interval version of the idempotent linear algebra and the corresponding 
universal algorithms are discussed). 

As a rule, iterative algorithms (beginning with the successive ap- 
proximation method) for solving differential equations (e.g., methods 
of Euler, Euler-Cauchy, Runge-Kutta, Adams, a number of important 
versions of the difference approximation method, and the like), meth- 
ods for calculating elementary and some special functions based on the 
expansion in Taylor's series and continuous fractions (Pade approxi- 
mations) and others are independent of the computer representation of 
numbers. 

3. Universal algorithms and accuracy of computations 

Calculations on computers usually are based on a floating-point 
arithmetic with a mantissa of a fixed length; i.e., computations are 
performed with fixed accuracy. Broadly speaking, with this approach 
only the relative rounding error is fixed, which can lead to a drastic loss 
of accuracy and invalid results (e.g., when summing series and subtract- 
ing close numbers). On the other hand, this approach provides rather 
high speed of computations. Many important numerical algorithms 
are designed to use floating-point arithmetic (with fixed accuracy) and 
ensure the maximum computation speed. However, these algorithms 
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are not universal. The above mentioned Gauss- Jacobi quadrature for- 
mulas, computation of elementary and special functions on the basis 
of the best polynomial or rational approximations or Pade-Chebyshev 
approximations, and some others belong to this type. Such algorithms 
use nontrivial constants specified with fixed accuracy. 

Recently, problems of accuracy, reliability, and authenticity of com- 
putations (including the effect of rounding errors) have gained much 
atention; in part, this fact is related to the ever-increasing performance 
of computer hardware. When errors in initial data and rounding errors 
strongly affect the computation results, such as in ill-posed problems, 
analysis of stability of solutions, etc., it is often useful to perform com- 
putations with improved and variable accuracy. In particular, the ratio- 
nal arithmetic, in which the rounding error is specified by the user [44], 
can be used for this purpose. This arithmetic is a useful complement to 
the interval analysis [54] . The corresponding computational algorithms 
must be universal (in the sense that they must be independent of the 
computer representation of numbers). 

4. Mathematics of semirings 

A broad class of universal algorithms is related to the concept of a 
semiring. We recall here the definition of a semiring (see, e.g., [21, 22]). 

4.1. Basic definitions. Consider a semiring, i.e., a set S endowed 
with two associative operations: addition © and multiplication such 
that addition is commutative, multiplication distributes over addition 
from either side, (resp., 1) is the neutral element of addition (resp., 
multiplication), x = x = for all x G S, and 0^1. Let 
the semiring S be partially ordered by a relation =$ such that is the 
least element and the inequality x =$ y implies that x © z =<; y © z, 
x z =4 V © z, and z x =4 z y for all x, y, z G S; in this case the 
semiring S is called positive (see, e.g., [22]). 

A semiring S is called a semifield if every nonzero element is invert- 
ible. 

A semiring S is called idempotent if x © x = x for all x G S, see, e.g., 
[4, 5, 7, 9, 10, 21-24, 31, 37-43, 48-53]. In this case the addition © 
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defines a canonical partial order ==! on the semiring S by the rule: 
x =^ y iff x © y = y. It is easy to prove that any idempotent semiring 
is positive with respect to this order. Note also that x © y — sup{a;, y} 
with respect to the canonical order. In the sequel, we shall assume 
that all idempotent semirings are ordered by the canonical partial order 
relation. 

We shall say that a positive (e.g., idempotent) semiring S is complete 
if it is complete as an ordered set. This means that for every subset 
T G S there exist elements sup T G S and inf T G S. 

The most well-known and important examples of positive semirings 
are "numerical" semirings consisting of (a subset of) real numbers and 
ordered by the usual linear order $C on R: the semiring R + with the 
usual operations © = +,© = • and neutral elements = 0, 1 = 1, the 
semiring R max = RU {—00} with the operations © = max, = + and 
neutral elements = —00, 1 = 0, the semiring R max = Rmax U {00}, 
where x =3! 00, x © 00 = 00 for all x, x 00 = 00 x = 00 if x ^ 0, 
and 00 = 00 0, and the semiring ^jnaxmin = [a, b], where —00 ^ 
a < b ^ +00, with the operations © = max, = min and neutral 
elements = a, 1 = b. The semirings R max , R max , and S [ ^ min = [a, b] 
are idempotent. The semirings R ma x, ^maimm; R+ = R+ U{°°} are 
complete. Remind that every partially ordered set can be imbedded to 
its completion (a minimal complete set containing the initial one). 

The semiring R m i n = R|J{oo} with operations © = min and = + 
and neutral elements = 00, 1 = is isomorphic to R ma x- 

The semiring R max is also called the max-plus algebra. The semifields 
Rmax and R m in are called tropical algebras. The term "tropical" initially 
appeared in [68] for a discrete version of the max-plus algebra as a 
suggestion of Ch. Choffrut, see also [26, 55, 61, 72]. 

Many mathematical constructions, notions, and results over the fields 
of real and complex numbers have nontrivial analogs over idempotent 
semirings. Idempotent semirings have become recently the object of 
investigation of new branches of mathematics, idempotent mathemat- 
ics and tropical geometry, see, e.g. [5, 13, 15-19, 35-28, 31, 32, 37-46, 
48-53, 55, 71, 72]. 
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4.2. Closure operations. Let a positive semiring S be endowed with 
a partial unary closure operation * such that x 4 y implies x* 4 V* 
and x* = 1 © (x* x) = 1 (B {x Q x*) on its domain of definition. 
In particular, 0* = 1 by definition. These axioms imply that x* = 
1 © x © x 2 © • • • © (x* a: ra ) if n ^ 1. Thus x* can be considered as a 
'regularized sum' of the series x* = 1 © x © x 2 © ... ; in an idempotent 
semiring, by definition, x* = sup{l, x, x 2 , . . . } if this supremum exists. 
So if S is complete, then the closure operation is well-defined for every 
element x G S. 

In numerical semirings the operation * is defined as follows: x* = 
(1 — x)^ 1 if x -< 1 in R + , or R + and x* = oo if rr )p 1 in R + ; x* = 1 
if x =<! 1 in R max and R max , x* = oo if x >- 1 in R max , a;* = 1 for all 
x in S^'-l „,.. In all other cases x* is undefined. Note that the closure 
operation is very easy to implement. 

4.3. Matrices over semirings. Denote by M&t mn (S) a set of all ma- 
trices A = (ctij) with m rows and n columns whose coefficients belong 
to a semiring S. The sum A © B of matrices A,B G Mat m „(5 r ) and the 
product AB of matrices A G Mat; m (S') and i? G Mat mn (S') are defined 
according to the usual rules of linear algebra: A © B = (a^- © fry) G 
Mat mn (S') and 

(m \ 

QaijObkA eM&t ln (S), 

where A G Mat; m (5 r ) and I? G Mat mn (S l ). Note that we write A5 
instead of A £>. 

If the semiring 5 is positive, then the set Mai mn (S) is ordered by 
the relation A = (oy) =$ B = (by) iff a^ =<! 6^- in 5 for all 1 ^ « ^ m, 
1 ^ j ^ n. 

The matrix multiplication is consistent with the order =<! in the fol- 
lowing sense: if A, A' G M&ti m (S), 5,5' G Mat mn (5 r ) and A ^ A', 
B 4 B', then AB 4 A'B' in Mat Zn (5 r ). The set Mat nn (5 r ) of square 
(n x n) matrices over a [positive, idempotent] semiring 5 forms a [posi- 
tive, idempotent] semiring with a zero element O = (0y), where Oij = 0, 
1 ^ hj ^ n, and a unit element E = (Sij), where 5ij = 1 if % = j and 
5ij = otherwise. 
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The set Mat nn is an example of a noncommutative semiring if n > 1. 

The closure operation in matrix semirings over a positive semiring S 
can be defined inductively (another way to do that see in [22] and 
below): A* = (an)* = (on) in Matn(S') and for any integer n > 1 and 
any matrix 

A = ( An Au 
\A 21 A 22/ 

where A u e M&t kk (S), A 12 G Mat fen _ fc (S), A 21 G Mat n _ fefe (S), A 22 G 
Mat n _fcn-fc(5'), 1 ^ k ^ n, by defintion, 

(1) A* = 

\ D*A 21 A* n D* 

where D = A 22 © A 21 Al ± Ai 2 . It can be proved that this definition of 
A* implies that the equality A* = A* A © E is satisfied and thus A* is 
a 'regularized sum' of the series E © A © A 2 © ... . 

Note that this recurrence relation coincides with the formulas of 
escalator method of matrix inversion in the traditional linear algebra 
over the field of real or complex numbers, up to the algebraic operations 
used. Hence this algorithm of matrix closure requires a polynomial 
number of operations in n. 

4.4. Discrete stationary Bellman equations. Let S be a positive 
semiring. The discrete stationary Bellman equation has the form 

(2) X = AX © B, 

where A e MaX nn (S), X, B G Mat ns (S), and the matrix X is unknown. 
Let A* be the closure of the matrix A. It follows from the identity 
A* = A* A © E that the matrix A*B satisfies this equation; moreover, 
it can be proved that for positive semirings this solution is the least in 
the set of solutions to equation (2) with respect to the partial order in 
Mat ns (S). 

4.5. Weighted directed graphs and matrices over semirings. 

Suppose that S is a semiring with zero and unity 1. It is well-known 
that any square matrix A = (a^) G M&t nn (S) specifies a weighted 
directed graph. This geometrical construction includes three kinds of 
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objects: the set X of n elements X\, . . . , x n called nodes, the set T of all 
ordered pairs (xi,Xj) such that a^ 7^ called arcs, and the mapping 
A: r — » S such that A(xi,Xj) = a„. The elements a^ of the semiring 
5 are called weights of the arcs. 

Conversely, any given weighted directed graph with n nodes specifies 
a unique matrix A G Mat nn (S l ). 

This definition allows for some pairs of nodes to be disconnected if 
the corresponding element of the matrix A is and for some channels 
to be "loops" with coincident ends if the matrix A has nonzero diag- 
onal elements. This concept is convenient for analysis of parallel and 
distributed computations and design of computing media and networks 
(see, e.g., [4, 53, 73]). 

Recall that a sequence of nodes of the form 

P = (yo,yi,---,Vk) 

with k ^ and (y^ yi+i) G T, i — 0, . . . , k — 1, is called a path of length 
k connecting y with y^. Denote the set of all such paths by Pk(yo, yu)- 
The weight A(p) of a path p G Pk{yo,yk) is defined to be the product 
of weights of arcs connecting consecutive nodes of the path: 

A(p) = A(y , Vl ) © • • • © A(j/ fc _i, y k ). 

By definition, for a 'path' p G Po(xi,Xj) of length k = the weight is 
1 if % — j and otherwise. 

For each matrix A G Mat nn (S') define A = E = (5{j) (where 5ij = 1 
Hi — j and 5ij = otherwise) and A k = AA k ^ 1 , k ^ 1. Let a\, be the 
(i, j)th element of the matrix A k . It is easily checked that 

Vl7 a Wl © ' ' ' © a ik-iik- 



a tJ 



i =i,i k =j 



Thus a]- is the supremum of the set of weights corresponding to all 
paths of length k connecting the node Xi = Xj with Xi k = Xj. 

Denote the elements of the matrix A* by a\* , i,j = l,...,n; then 



"1:1 



<>■■■■= © © MP)- 



,(*) _ 

0^fc<oo p&Pk(xi,Xj) 
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The closure matrix A* solves the well-known algebraic path prob- 
lem, which is formulated as follows: for each pair (xi,Xj) calculate the 
supremum of weights of all paths (of arbitrary length) connecting node 
Xi with node Xj. The closure operation in matrix semirings has been 
studied extensively (see, e.g., [1, 2, 5-7, 10, 11, 15-17, 22-26, 31, 32, 
46] and references therein). 

Example 1. The shortest path problem. Let S = R m i n , so the 
weights are real numbers. In this case 

A(p) = A(y , y x ) + A(y 1: y 2 ) H h A(y k _ 1: y k ). 

If the element a^ specifies the length of the arc (xi,Xj) in some metric, 
then a)j is the length of the shortest path connecting Xi with Xj. 
Example 2. The maximal path width problem. Let S = R U 

{0, 1} with © = max, = min. Then 

aff = max A(p), A(p) = min(A(y , j/i), . . . , A(y k _ 1} y k )). 
pe U p k(xi,xj) 

If the element a^- specifies the "width" of the arc (xi,Xj), then the 
width of a path p is defined as the minimal width of its constituting 
arcs and the element a\ • gives the supremum of possible widths of all 
paths connecting Xi with Xj. 

Example 3. A simple dynamic programming problem. Let S = 
R-max and suppose a„ gives the profit corresponding to the transition 
from Xi to Xj. Define the vector B = (6j) G Mat„i(R max ) whose element 
bi gives the terminal profit corresponding to exiting from the graph 
through the node X{. Of course, negative profits (or, rather, losses) 
are allowed. Let m be the total profit corresponding to a path p € 

m = A(p) + by 

Then it is easy to check that the supremum of profits that can be 
achieved on paths of length k beginning at the node Xi is equal to 
(A k B)i and the supremum of profits achievable without a restriction 
on the length of a path equals (A*B)i. 

Example 4. The matrix inversion problem. Note that in the for- 
mulas of this section we are using distributivity of the multiplication 
with respect to the addition © but do not use the idempotency axiom. 
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Thus the algebraic path problem can be posed for a nonidempotent 
semiring S as well (see, e.g., [66]). For instance, if S = R, then 

A* = E + A + A 2 + ■ ■ ■ = (E - A)- 1 . 

If || A || > 1 but the matrix E — A is invertible, then this expression 
defines a regularized sum of the divergent matrix power series ^2 i>0 A 1 . 

There are many other important examples of problems (in different 
areas) related to algorithms of linear algebra over semirings (transitive 
closures of relations, accessible sets, critical paths, paths of greatest 
capacities, the most reliable paths, interval and other problems), see 
[1, 2, 4, 5, 10-13, 15-18, 20, 22-27, 31, 32, 45, 46, 53, 57, 58, 63-68, 
74-77]. 

We emphasize that this connection between the matrix closure op- 
eration and solution to the Bellman equation gives rise to a number 
of different algorithms for numerical calculation of the closure matrix. 
All these algorithms are adaptations of the well-known algorithms of 
the traditional computational linear algebra, such as the Gauss-Jordan 
elimination, various iterative and escalator schemes, etc. This is a spe- 
cial case of the idempotent superposition principle (see below). 

In fact, the theory of the discrete stationary Bellman equation can 
be developed using the identity A* = AA* © E as an additional axiom 
without any substantive interpretation (the so-called closed semirings, 
see, e.g., [7, 22, 36, 66]). 

5. Universal algorithms of linear algebra over semirings 

The most important linear algebra problem is to solve the system of 
linear equations 

(3) AX = B, 

where A is a matrix with elements from the basic field and X and B are 
vectors (or matrices) with elements from the same field. It is required 
to find X if A and B are given. If A in (3) is not the identity matrix 
/, then system (3) can be written in form (2), i.e., 

X = AX + B. (2') 
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It is well known that the form (2) or (2') is convenient for using the suc- 
cessive approximation method. Applying this method with the initial 
approximation X = 0, we obtain the solution 

(4) X = A*B, 
where 

(5) A* = I + A + A 2 + • • • + A n + • • • 
On the other hand, it is clear that 

(6) A* = (I-A)-\ 

if the matrix I — A is invertible. The inverse matrix (I — A)~ x can be 
considered as a regularized sum of the formal series (5). 

The above considerations can be extended to a broad class of semir- 
ings. 

The closure operation for matrix semirings Mat n (S') can be defined 
and computed in terms of the closure operation for S (see section 4.3 
above); some methods are described in [1, 2, 7, 10, 11, 22, 23-25, 
32, 35, 42, 46, 65-67]. One such method is described below (LDM- 
factorization). 

If S is a field, then, by definition, x* = (1 — x)~ l for any x ^ 1. If S 
is an idempotent semiring, then, by definition, 

(7) x* = 1 © x © x 2 © • • • = sup{l, x, x 2 ,... }, 

if this supremum exists. Recall that it exists if S is complete, see 
section 4.2. 

Consider a nontrivial universal algorithm applicable to matrices over 
semirings with the closure operation defined. 

Example 5: Semiring LDM -Factorization. 

Factorization of a matrix into the product A = LDM, where L 
and M are lower and upper triangular matrices with a unit diagonal, 
respectively, and D is a diagonal matrix, is used for solving matrix 
equations AX = B. We construct a similar decomposition for the 
Bellman equation X = AX © B. 
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For the case AX = B, the decomposition A = LDM induces the 
following decomposition of the initial equation: 

(8) LZ = B, DY = Z, MX = Y. 
Hence, we have 

(9) A' 1 = M^D^L- 1 , 

if A is invertible. In essence, it is sufficient to find the matrices L, D 
and M, since the linear system (8) is easily solved by a combination 
of the forward substitution for Z, the trivial inversion of a diagonal 
matrix for Y, and the back substitution for X. 
Using (8) as a pattern, we can write 

(10) Z = LZ®B, Y = DY®Z, X = MX ®Y. 
Then 

(11) A* = M*D*L*. 

A triple (L, D, M) consisting of a lower triangular, diagonal, and upper 
triangular matrices is called an LD M- factorization of a matrix A if 
relations (10) and (11) are satisfied. We note that in this case, the 
principal diagonals of L and M are zero. 

The modification of the notion of LDM-factorization used in matrix 
analysis for the equation AX = B is constructed in analogy with the 
construct suggested by Carre in [10, 11] for L£/-factorization. 

We stress that the algorithm described below can be applied to ma- 
trix computations over any semiring under the condition that the unary 
operation a t— > a* is applicable every time it is encountered in the com- 
putational process. Indeed, when constructing the algorithm, we use 
only the basic semiring operations of addition © and multiplication 
and the properties of associativity, commutativity of addition, and 
distributivity of multiplication over addition. 

If A is a symmetric matrix over a semiring with a commutative mul- 
tiplication, the amount of computations can be halved, since M and L 
are mapped into each other under transposition. 

We begin with the case of a triangular matrix A = L (or A = M). 
Then, finding X is reduced to the forward (or back) substitution. 



MATHEMATICS OF SEMIRINGS 17 

Forward substitution 
We are given: 

• L = |U'-||r o -_i, where l l - = for i < j (a lower triangular matrix 
with a zero diagonal); 

• £= 116*11?=!. 

It is required to find the solution X = ||x J ||™ =1 to the equation X = 
LX © B. The program fragment solving this problem is as follows, 
for % = 1 to n do 
{ x* := b' 1 ; 

for j = 1 to % — 1 do 

x l := x* @ (l) O tf); } 

Back substitution 
We are given 

• M — ||m*||™ =1 , where rrij = for % > j (an upper triangular 
matrix with a zero diagonal); 

• B= 116*11?=!. 

It is required to find the solution X = ||a; l ||™ =1 to the equation X = 
MX © B. The program fragment solving this problem is as follows. 

for % — n to 1 step —1 do 
{ x* := 6*; 

for j — n to % + 1 step —1 do 

x 4 := a; 4 © (m* x*); } 

Both algorithms require (n 2 — n)/2 operations © and 0. 

Closure of a diagonal matrix 
We are given 

• D =diag(di,...,d n ); 

• B = 116*11?=!. 

It is required to find the solution X = ||x l ||™ =1 to the equation X = 
DX © B. The program fragment solving this problem is as follows. 

for % — 1 to n do 
x* := (dO* 6*'; 
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This algorithm requires n operations * and n multiplications 0. 

General case 
We are given 

• L = ||/}||" J=1 , where /*■ = if % < j; 

• D = di&g(d 1 ,...,d n ); 

• M = ||m*-||r,--i, where m\ = if i > ?': 

• B = 116*11?=!. 

It is required to find the solution X = ||x l ||™ =1 to the equation X = 
AX © B, where L, D, and M form the LDM-factorization of A. The 
program fragment solving this problem is as follows. 

FORWARD SUBSTITUTION 
for % — 1 to n do 

{ x l := b l ; 

for j — 1 to i — 1 do 

x l := x l @ (I) Q x j ); } 
CLOSURE OF A DIAGONAL MATRIX 
for i — 1 to n do 

x l ■= (di)* V] 
BACK SUBSTITUTION 
for % — n to 1 step —1 do 
{ for j — n to % + 1 step —1 do 

af := x* © (to*- x J '); } 

Note that x % is not initialized in the course of the back substitution. 
The algorithm requires n 2 — n operations ©, n 2 operations 0, and n 
operations *. 

LDM- factorization 

We are given 

_ 4 _ \\ n i \\n 
•A— 110^-11^-=!. 

It is required to find the LDM-factorization of A: L = ||/!||™, =1 , 
D = diag(di, . . . , d n ), and M = Hto*-!!™^, where /* = if i < j, and 
rrij = if i > j . 
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The program uses the following internal variables: 

_ n — \\ r i \\n 

• v = IMIJU; 

• d. 

INITIALISATION 

for % = 1 to n do 

for j = 1 to n do 

c i = a i' 



MAIN LOOP 

for j = 1 to n do 
{ for i = 1 to j do 



v l := ay, 



for k = 1 to j — 1 do 

for 2 = k + 1 to j do 

^ :=««©(o|0i; k ); 
for % — 1 to j — 1 do 





a* :=(aj)*0u»; 






a J : 


= ^"; 






for /c = 1 to j — 1 do 








for i = j + 1 to n 


do 






a} := aj © (a] 


i® 


w fc ); 


d = 


= (v j Y; 






for 


i = j + 1 to n do 







a i := a i ® ^' } 

This algorithm requires (2n 3 — 3n 2 +n)/6 operations ©, (2n 3 + 3n 2 — 
5n)/6 operations 0, and n(n + 1)/2 operations *. After its completion, 
the matrices L, D, and M are contained, respectively, in the lower 
triangle, on the diagonal, and in the upper triangle of the matrix C . 
In the case when A is symmetric about the principal diagonal and the 
semiring over which the matrix is defined is commutative, the algorithm 
can be modified in such a way that the number of operations is reduced 
approximately by a factor of two. 

Other examples can be found in [10, 11, 22-25, 35, 36, 66, 67]. 
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Note that to compute the matrices A* and A*B it is convenient to 
solve the Bellman equation (2). 

Some other interesting and important problems of linear algebra over 
semirings are examined, e.g., in [8, 9, 12, 18, 20, 22-25, 27, 57, 58, 59, 

60, 74-77]. 

6. The idempotent correspondence principle 

There is a nontrivial analogy between mathematics of semirings 
and quantum mechanics. For example, the field of real numbers can 
be treated as a "quantum object" with respect to idempotent semir- 
ings. So idempotent semirings can be treated as "classical" or "semi- 
classical" objects with respect to the field of real numbers. 

Let R be the field of real numbers and R + the subset of all non- 
negative numbers. Consider the following change of variables: 

u i— > w — h In u, 

where u G R+ \ {0}, h > 0; thus u = e w ^ h , w G R. Denote by the 
additional element — oo and by S the extended real line R U {0}. The 
above change of variables has a natural extension D^ to the whole S 
by Dh(0) = 0; also, we denote -D^(l) =0 = 1. 

Denote by Sh the set S equipped with the two operations Q)h (gen- 
eralized addition) and 0^ (generalized multiplication) such that Dh 
is a homomorphism of {R + ,+,-} to {S, ®h, &h}- This means that 
D h (ui + u 2 ) = D h (m) ® h D h (u 2 ) and D h (u x ■ u 2 ) = D h (m) & h D h (u 2 ), 
i.e., wi &h w 2 = w\ + w 2 and w\ (Bh w 2 = hln(e Wl ^ h + e W2 ^ h ). It is easy 
to prove that w 1 Q)h u) 2 — > max{wi, w 2 } as h — > 0. 

Denote by R max the set S = R U {0} equipped with operations 
© = max and = +, where = — oo,l = 0as above. Algebraic 
structures in R + and Sh are isomorphic; therefore R max is a result of 
a deformation of the structure in R + . 

We stress the obvious analogy with the quantization procedure, where 
h is the analog of the Planck constant. In these terms, R + (or R) 
plays the part of a "quantum object" while R max acts as a "classical" 
or "semi-classical" object that arises as the result of a dequantization 
of this quantum object. 
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Likewise, denote by R m i n the set R U {0} equipped with operations 
© = min and = +, where = +oo and 1 = 0. Clearly, the 
corresponding dequantization procedure is generated by the change of 
variables «4w= — hlnu. 

Consider also the set RU {0, 1}, where = — oo, 1 = +oo, together 
with the operations © = max and = min. Obviously, it can be 
obtained as a result of a "second dequantization" of R or R + . 

There is a natural transition from the field of real numbers or com- 
plex numbers to the semiring R max (or R mm ). This is a composition of 
the mapping x *-¥ \x\ and the deformation described above. 

In general an idempotent dequantization is a transition from a basic 
field to an idempotent semiring in mathematical concepts, construc- 
tions and results, see [38-40] for details. 

For example, the basic object of the traditional calculus is a function 
defined on some set X and taking its values in the field R (or C); 
its idempotent analog is a map X — >■ S, where X is some set and 
S = Rmin, Rmax, or another idempotent semiring. Let us show that 
redefinition of basic constructions of traditional calculus in terms of 
idempotent mathematics can yield interesting and nontrivial results 
(see, e.g., [32, 37-43, 45, 46, 48-53, 71, 72], for details, additional 
examples and generalizations). 

Example 6. Integration and measures. To define an idempotent 
analog of the Riemann integral, consider a Riemann sum for a function 
(p(x), x G X = [a,b], and substitute semiring operations © and for 
operations + and • (usual addition and multiplication) in its expression 
(for the sake of being definite, consider the semiring R ma x): 

N N 

y^(p(xi)-Ai h-> ffi^)0Ai= max (y9(x;) + A;), 

^— • ^-^ i=l,...,N 

i=l i=l 

where a = x < X\ < • • • < x N = b, Aj = x, — Xj_i, i — 1, . . . , N. As 
maxj Aj — ¥ 0, the integral sum tends to 



<p(x) dx = sup <p(x) 

X x&X 
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for any function ip: X — > R max that is bounded. In general, for any set 
X the set function 

m<p(B) = sup(p(x), B C X, 
xeB 

is called an R max -meas«re on X; since rn ip ([j a B a ) = sup a m,p(B a ), this 
measure is completely additive. An idempotent integral with respect 
to this measure is defined as 



ip(x) drrip = / ip(x) ip(x) dx = sup (ip(x) + ip(x)). 
x ' Jx xex 

Using the standard partial order it is possible to generalize these 
definitions for the case of arbitrary idempotent semirings. 
Example 7. Fourier Legendre transform. Consider the topolog- 
ical group G = R n . The usual Fourier-Laplace transform is defined 

as 

<p(x) ^ £(£) = / e*' x <p(x) dx, 
Jg 

where exp(i£ • x) is a character of the group G, i.e., a solution of the 
following functional equation: 

f(x + y) = f(x)f(y). 

The idempotent analog of this equation is 

f(x + y) = f(x)&f(y) = f(x) + f(y). 

Hence natural "idempotent characters" of the group G are linear func- 
tions of the form x i— ¥ £-x = £1^1 + ' ' ' J r£,n x n- Thus the Fourier-Laplace 
transform turns into 



<p(x) i->- </?(£) — £ ■ x <p(x) dx = sup (£ • x + <p(x)). 
Jg xeG 

This is the well-known Legendre (or Fenchel) transform. Examples re- 
lated to an important version of matrix algebra are discussed in section 
4 above. 

These examples suggest the following formulation of the idempotent 
correspondence principle [39, 40]: 

There exists a heuristic correspondence between interest- 
ing, useful, and important constructions and results over 
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Figure 1. Relations between idempotent and tradi- 
tional mathematics. 

the field of real (or complex) numbers and similar con- 
structions and results over idempotent semirings in the 
spirit of N. Bohr's correspondence principle in quantum 
mechanics. 

So idempotent mathematics can be treated as a "classical shadow 
(or counterpart)" of the traditional Mathematics over fields. 

A systematic application of this correspondence principle leads to a 
variety of theoretical and applied results, see, e.g. [37-43, 45, 46, 55, 
71, 72]. Relations between idempotent and traditional mathematics 
are presented in Fig 1. Relations to quantum physics are discussed in 
detail, e.g., in [38]. 



7. The superposition principle and parallel computing 

In quantum mechanics the superposition principle means that the 
Schrodinger equation (which is basic for the theory) is linear. Simi- 
larly in idempotent mathematics the idempotent superposition prin- 
ciple means that some important and basic problems and equations 
(e.g., optimization problems, the Bellman equation and its versions 
and generalizations, the Hamilton-Jacobi equation) nonlinear in the 
usual sense can be treated as linear over appropriate idempotent semir- 
ings (in a general form this superposition principle was formulated by 
V. P. Maslov), see [48-53]. 
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Linearity of the Hamilton- Jacobi equation over R m i n (and R max ) 
can be deduced from the usual linearity (over C) of the correspond- 
ing Schrodinger equation by means of the dequantization procedure 
described above (in Section 4). In this case the parameter h of this 
dequantization coincides with ifi , where h~ is the Planck constant; so 
in this case h must take imaginary values (because h > 0; see [38] for 
details). Of course, this is closely related to variational principles of 
mechanics. 

The situation is similar for the differential Bellman equation (see, 
e.g., [32]) and the discrete version of the Bellman equations, see section 
4 above. 

It is well known that linear problems and equations are especially 
convenient for parallelization, see, e.g., [73]. Standard methods (in- 
cluding the so-called block methods) constructed in the framework of 
the traditional mathematics can be extended to universal algorithms 
over semirings (the correspondence principle!). For example, formula 
(1) discussed in section 4.3 leads to a simple block method for paral- 
lelization of the closure operations. Other standard methods of linear 
algebra [73] can be used in a similar way. 

8. The correspondence principle for computations 

Of course, the idempotent correspondence principle is valid for algo- 
rithms as well as for their software and hardware implementations [39, 
40, 42, 43]. Thus: 

If we have an important and interesting numerical al- 
gorithm, then there is a good chance that its semiring 
analogs are important and interesting as well. 

In particular, according to the superposition principle, analogs of 
linear algebra algorithms are especially important. Note that numer- 
ical algorithms for standard infinite-dimensional linear problems over 
idempotent semirings (i.e., for problems related to idempotent inte- 
gration, integral operators and transformations, the Hamilton- Jacobi 
and generalized Bellman equations) deal with the corresponding finite- 
dimensional (or finite) "linear approximations". Nonlinear algorithms 
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often can be approximated by linear ones. Thus the idempotent linear 
algebra is a basis for the idempotent numerical analysis. 

Moreover, it is well-known that linear algebra algorithms easily lend 
themselves to parallel computation; their idempotent analogs admit 
parallelization as well. Thus we obtain a systematic way of applying 
parallel computing to optimization problems. 

Basic algorithms of linear algebra (such as inner product of two vec- 
tors, matrix addition and multiplication, etc.) often do not depend on 
concrete semirings, as well as on the nature of domains containing the 
elements of vectors and matrices. Algorithms to construct the closure 
A* = l®A®A 2 ®---®A n ®--- = 0^° =1 A n of an idempotent ma- 
trix A can be derived from standard methods for calculating (1 — A)^ 1 . 
For the Gauss- Jordan elimination method (via LU-decomposition) this 
trick was used in [66], and the corresponding algorithm is universal and 
can be applied both to the Bellman equation and to computing the in- 
verse of a real (or complex) matrix (1 — A). Computation of A~ l can 
be derived from this universal algorithm with some obvious cosmetic 
transformations. 

Thus it seems reasonable to develop universal algorithms that can 
deal equally well with initial data of different domains sharing the same 
basic structure [39, 40, 43]. 



9. The correspondence principle for hardware design 

A systematic application of the correspondence principle to computer 
calculations leads to a unifying approach to software and hardware 
design. 

The most important and standard numerical algorithms have many 
hardware realizations in the form of technical devices or special proces- 
sors. These devices often can be used as prototypes for new hardware 
units generated by substitution of the usual arithmetic operations for 
its semiring analogs and by addition tools for performing neutral ele- 
ments and 1 (the latter usually is not difficult). Of course, the case of 
numerical semirings consisting of real numbers (maybe except neutral 
elements) and semirings of numerical intervals is the most simple and 
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natural [38-43, 45, 46, 69]. Note that for semiflelds (including R max 
and R m i n ) the operation of division is also defined. 

Good and efficient technical ideas and decisions can be transposed 
from prototypes into new hardware units. Thus the correspondence 
principle generated a regular heuristic method for hardware design. 
Note that to get a patent it is necessary to present the so-called 'inven- 
tion formula', that is to indicate a prototype for the suggested device 
and the difference between these devices. 

Consider (as a typical example) the most popular and important 
algorithm of computing the scalar product of two vectors: 

(12) (X, y) = XiJ/i + x 2 y 2 H Vx n y n . 

The universal version of (12) for any semiring A is obvious: 

(13) (x, y) = (xi yi) © (x 2 y 2 ) © • • • © (x n & y n ). 

In the case A = R max this formula turns into the following one: 

(14) (x, y) = max{xi + y u x 2 + y 2 , • • • , x n + y n }. 

This calculation is standard for many optimization algorithms, so it 
is useful to construct a hardware unit for computing (14). There are 
many different devices (and patents) for computing (12) and every such 
device can be used as a prototype to construct a new device for comput- 
ing (14) and even (13). Many processors for matrix multiplication and 
for other algorithms of linear algebra are based on computing scalar 
products and on the corresponding "elementary" devices respectively, 
etc. 

There are some methods to make these new devices more universal 
than their prototypes. There is a modest collection of possible opera- 
tions for standard numerical semirings: max, min, and the usual arith- 
metic operations. So, it is easy to construct programmable hardware 
processors with variable basic operations. Using modern technologies 
it is possible to construct cheap special-purpose multi-processor chips 
implementing examined algorithms. The so-called systolic processors 
are especially convenient for this purpose. A systolic array is a 'homo- 
geneous' computing medium consisting of elementary processors, where 
the general scheme and processor connections are simple and regular. 
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Every elementary processor pumps data in and out performing elemen- 
tary operations in a such way that the corresponding data flow is kept 
up in the computing medium; there is an analogy with the blood cir- 
culation and this is a reason for the term "systolic", see e.g., [39, 40, 
43, 52, 65-67]. 

Concrete systolic processors for the general algebraic path problem 
are presented in [65-67]. In particular, there is a systolic array of 
n(n + 1) elementary processors which performs computations of the 
Gauss-Jordan elimination algorithm and can solve the algebraic path 
problem within 5n— 2 time steps. Of course, hardware implementations 
for important and popular basic algorithms increase the speed of data 
processing. 

The so-called GPGPU (General-Purpose computing on Graphics 
Processing Units) technique is another important field for applications 
of the correspondence principle. The thing is that graphic processing 
units (hidden in modern laptop and desktop computers) are potentially 
powerful processors for solving numerical problems. The recent tremen- 
dous progress in graphical processing hardware and software leads to 
new "open" programmable parallel computational devices (special pro- 
cessors), see, e.g., [78-80]. These devices are going to be standard for 
coming PC (personal computers) generations. Initially used for graph- 
ical processing only (at that time they were called GPU), today they 
are used for various fields, including audio and video processing, com- 
puter simulation, and encryption. But this list can be considerably 
enlarged following the correspondence principle: the basic operations 
would be used as parameters. Using the technique described in this 
paper (see also our references), standard linear algebra algorithms can 
be used for solving different problems in different areas. In fact, the 
hardware supports all operations needed for the most important idem- 
potent semirings: plus, times, min, max. The most popular linear alge- 
bra packages [ATLAS (Automatically Tuned Linear Algebra Software), 
LAPACK, PLASMA (Parallel Linear Algebra for Scalable Multicore 
Architectures)] can already use GPGPU, see [81-83]. We propose to 
make these tools more powerful by using parameterized algorithms. 
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Linear algebra over the most important numerical semirings gener- 
ates solutions for many concrete problems in different areas, see, e.g., 
sections 4.4 and 4.5 above and references indicated in these sections. 

Note that to be consistent with operations we have to redefine zero 
(0) and unit (1) elements (see above); comparison operations must 
be also redefined as it is described in section 4.1 "Basic definitions." 
Once the operations are redefined, then the most of basic linear algebra 
algorithms, including back and forward substitution, Gauss elimination 
method, Jordan elimination method and others could be rewritten for 
new domains and data structures. Combined with the power of the new 
parallel hardware this approach could change PC from entertainment 
devices to power full instruments. 

10. The correspondence principle for software design 

Software implementations for universal semiring algorithms are not 
as efficient as hardware ones (with respect to the computation speed) 
but they are much more flexible. Program modules can deal with ab- 
stract (and variable) operations and data types. Concrete values for 
these operations and data types can be defined by the corresponding 
input data. In this case concrete operations and data types are gener- 
ated by means of additional program modules. For programs written 
in this manner it is convenient to use special techniques of the so- 
called object oriented (and functional) design, see, e.g., [47, 62, 70]. 
Fortunately, powerful tools supporting the object-oriented software de- 
sign have recently appeared including compilers for real and convenient 
programming languages (e.g. C ++ and Java) and modern computer al- 
gebra systems. 

Recently, this type of programming technique has been dubbed generic 
programming (see, e.g., [6, 62]). To help automate the generic program- 
ming, the so-called Standard Template Library (STL) was developed in 
the framework of C ++ [62, 70]. However, high-level tools, such as STL, 
possess both obvious advantages and some disadvantages and must be 
used with caution. 

It seems that it is natural to obtain an implementation of the corre- 
spondence principle approach to scientific calculations in the form of a 
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powerful software system based on a collection of universal algorithms. 
This approach ensures a working time reduction for programmers and 
users because of the software unification. The arbitrary necessary ac- 
curacy and safety of numeric calculations can be ensured as well. 

The system has to contain several levels (including programmer and 
user levels) and many modules. 

Roughly speaking, it must be divided into three parts. The first part 
contains modules that implement domain modules (finite representa- 
tions of basic mathematical objects). The second part implements uni- 
versal (invariant) calculation methods. The third part contains mod- 
ules implementing model dependent algorithms. These modules may 
be used in user programs written in C ++ , Java, Maple, Mathlab etc. 

The system has to contain the following modules: 

— Domain modules: 

— infinite precision integers; 

— rational numbers; 

— finite precision rational numbers (see [44]); 

— finite precision complex rational numbers; 

— fixed- and floating-slash rational numbers; 

— complex rational numbers; 

— arbitrary precision floating-point real numbers; 

— arbitrary precision complex numbers; 

— p-adic numbers; 

— interval numbers; 

— ring of polynomials over different rings; 

— idempotent semirings; 

— interval idempotent semirings; 

— and others. 

— Algorithms: 

— linear algebra; 

— numerical integration; 

— roots of polynomials; 

— spline interpolations and approximations; 

— rational and polynomial interpolations and approximations; 
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— special functions calculation; 

— differential equations; 

— optimization and optimal control; 

— idempotent functional analysis; 

— and others. 

This software system may be especially useful for designers of algo- 
rithms, software engineers, students and mathematicians. 

Note that there are some software systems oriented to calculations 
with idempotent semirings like R max ; see, e.g., [64]. However these 
systems do not support universal algorithms. 

11. Interval analysis in idempotent mathematics 

Traditional interval analysis is a nontrivial and popular mathemat- 
ical area, see, e.g., [3, 20, 34, 54, 56, 59]. An "idempotent" version 
of interval analysis (and moreover interval analysis over positive semir- 
ings) appeared in [45, 46, 69]. Later appeared rather many publications 
on the subject, see, e.g., [12, 20, 27, 57, 58, 77]. Interval analysis over 
the positive semiring R + was discussed in [8]. 

Let a set S be partially ordered by a relation =<;. A closed interval 
in S is a subset of the form x= [x, x] = {xG5'| x^x^x}, where 
the elements x =<! x are called lower and upper bounds of the interval x. 
The order =^ induces a partial ordering on the set of all closed intervals 
in S: x =<! y iff x =<! y and x =<! y. 

A weak interval extension I(S) of a positive semiring S is the set 
of all closed intervals in S endowed with operations © and defined 
as x © y = [x © y, x © y], x y = [x0y,x0y] and a partial order 
induced by the order in S. The closure operation in I(S) is defined by 
x* = [x*,x*]. There are some other interval extensions (including the 
so-called strong interval extension [46] ) but the weak extension is more 
convenient. 

The extension I(S) is positive; I(S) is idempotent if S is an idempo- 
tent semiring. A universal algorithm over S can be applied to I(S) and 
we shall get an interval version of the initial algorithm. Usually both 
the versions have the same complexity. For the discrete stationary Bell- 
man equation and the corresponding optimization problems on graphs 
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interval analysis was examined in [45, 46] in details. Other problems 
of idempotent linear algebra were examined in [12, 27, 57, 58, 77]. 

Idempotent mathematics appears to be remarkably simpler than its 
traditional analog. For example, in traditional interval arithmetic, mul- 
tiplication of intervals is not distributive with respect to addition of in- 
tervals, whereas in idempotent interval arithmetic this distributivity is 
preserved. Moreover, in traditional interval analysis the set of all square 
interval matrices of a given order does not form even a semigroup with 
respect to matrix multiplication: this operation is not associative since 
distributivity is lost in the traditional interval arithmetic. On the con- 
trary, in the idempotent (and positive) case associativity is preserved. 
Finally, in traditional interval analysis some problems of linear algebra, 
such as solution of a linear system of interval equations, can be very 
difficult (generally speaking, they are iVP-hard, see [14, 20, 33, 34] and 
references therein). It was noticed in [45, 46] that in the idempotent 
case solving an interval linear system requires a polynomial number of 
operations (similarly to the usual Gauss elimination algorithm). Two 
properties that make the idempotent interval arithmetic so simple are 
monotonicity of arithmetic operations and positivity of all elements of 
an idempotent semiring. 

Usually interval estimates in idempotent mathematics are exact. In 
the traditional theory such estimates tend to be overly pessimistic. 
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